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Abstract 

It has been recently proposed that the matter-antimatter asymmetry of the universe 
may have its origin in "post-sphaleron baryogenesis" (PSB). It is a TeV scale mech- 
anism that is testable at the LHC and other low energy experiments. In this paper 
we present a theory of PSB within a quark-lepton unified scheme based on the gauge 
group SU(2)l x SU(2)r x SU(4) c that allows a direct connection between the baryon 
asymmetry and neutrino mass matrix. The flavor changing neutral current constraints 
on the model allow successful baryogenesis only for an inverted mass hierarchy for neu- 
trinos, which can be tested in the proposed long base line neutrino experiments. The 
model also predicts observable neutron-antineutron oscillation accessible to the next 
generation of experiments as well as TeV scale colored scalars within reach of LHC. 



1 Introduction 



Since Sakharov first suggested the three conditions that would have to be satisfied by a 
microphysical theory to generate matter-antimatter asymmetry of the universe p], many 
beyond the standard model scenarios have been constructed for this purpose. The very ear- 
liest ones that used proton decay in grand unified theories for this purpose run into difficulty 
on several counts: first is that successful inflation scenarios generally have reheating temper- 
atures which are below the generic baryogenesis temperatures, especially in the context of 
super symmetry, so that any GUT generated baryon number is erased by inflation; secondly, 
if baryogenesis is caused by B — L conserving interactions as in SU (5) models, they will be 
destroyed by electroweak sphalerons that are in equilibrium down to about 100 GeV. 

In the mid 80's, a new mechanism was suggested that uses baryogenesis via leptogene- 
sis [2]. This mechanism is very attractive since it arises within the framework of the seesaw 
mechanism [3j that explains small neutrino masses. Here the initial lepton asymmetry is 
created far below the GUT scale and is then converted by the electroweak sphalerons j?J to a 
baryon asymmetry. This mechanism depends crucially on the properties of the electroweak 
sphaleron [I] which serves as the source of B violation. While this is one of the most widely 
discussed schemes in literature today [5], it may also have problems since adequate leptoge- 
nesis in these models implies a lower bound on the leptogenesis scale [6] which is above the 
allowed reheating scale in supersymmetric models [?]. It is also not so easy to test by low 
energy experiments. 

It is therefore important to explore alternative mechanisms that can explain the matter- 
antimatter asymmetry from particle decays around 100 GeV temperatures which do not 
conflict with the above bounds on reheat temperatures and at the same time yield testable 
consequences at LHC and other low energy experiments. Such a mechanism was proposed in 
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two recent papers [8j [9] , where it was shown that with the use of higher dimensional baryon 
violating operators, baryogenesis can occur after the electro-weak sphalerons have gone out 
of thermal equilibrium. This mechanism was called Post-Sphaleron Baryogenesis (PSB). 

In this paper, we propose a theory for this mechanism within a quark-lepton unified gauge 
model for neutrino masses based on the gauge group SU(2) l x SU (2) r x SU (4) c [10] using the 
symmetry breaking setup discussed in Ref. [II] . We show that in this model, quark-lepton 
unification allows us to relate the baryon asymmetry directly to neutrino masses via the 
type II seesaw mechanism [12J. The main result of our work is that successful baryogenesis 
can occur only for an inverted mass hierarchy for neutrinos with a relatively large #13. Both 
of these predictions will be tested in ongoing experiments searching for neutrino-less double 
beta decay and neutrino oscillations. 

The salient feature of PSB mechanism is that baryogenesis occurs via the direct decay 
of a scalar boson S r having a weak scale mass and a higher dimensional baryon violating 
coupling. S r is the real part of a baryon number carrying complex scalar S, which acquires 
a vacuum expectation value (vev). In the context of the SU{2) L x SU{2) R x SU(4) C model, 
S r is the real part of a Higgs scalar field belonging to (1, 3, 10) representation of whose vev 
breaks the SU(2) R x SU(4) C symmetry down to the U(l)y of the standard model. The 
decays S r — > 6q c and S r — *> 6q c provide the source for B asymmetry. When the S'-field has a 
vev, the decay process generates an interaction that causes neutron-antineutron oscillation as 
shown in Ref. pj]. The parameter domain of our theory where adequate baryogenesis occurs 
predicts that neutron-antineutron oscillation should occur at a rate observable in currently 
available reactor facilities. 



3 



2 Basic ingredients of Post-Sphaleron Baryogenesis and 

its SU(2) L x SU(2) R x SU{4) C embedding 



A starting Lagrangian for PSB that gives rise to the higher dimensional 5-violating decay 
is given by [8] 

+—SA u c u cA d c d cA d c d c + —SA d c d cA^c d c + h.c. (1) 

Here the A u c u c, etc. are color sextet scalar fields. From the above equation, we see that 
when the scalar field S, which has B — L = 2 is given a vev, it leads to cubic scalar field 
couplings of the type A^cA^cA^^ and A d c dc A^ cdc leading to baryon number violation by 
two units. 

We note that not all of the (A u c u c, A u c d c, A d c d c) fields are needed for B violation and 
n <r-> n oscillation: either (A u c d c, A d c d c) or (A u o u c, A^c) pair will do. In fact, consistency 
with flavor changing neutral current constraints and n — n oscillation limits allow for only 
two of these three scalar states to be light near the TeV scale. The third state (in our case 
A u c u c, as we will see below) will have mass of order 100 TeV. 

Baryon asymmetry arises in this scheme from WMoop corrections to the S r decays and 
is therefore directly linked to CKM mixing [8]. 

The constraints on the parameter space of the model arise from the fact that the decay 
of S r occurs below 100 GeV and above 200 MeV or so - the former to ensure that the 
sphalerons do not play any role in baryogenesis and the latter so that quarks in the cosmic 
soup have not combined to form hadrons, which will affect the decay estimates - and from 
the fact that the model must reproduce observed neutrino masses and mixings. If baryon 
asymmetry is created above the electroweak phase transition temperature, all of the baryon 
asymmetry will be washed out since there are both B + L violating sphaleron interactions as 
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well as Wr mediated AL = 2 scatterings of right-handed Majorana neutrinos in equilibrium 
at that temperature. 

Before discussing the constraints on the parameters of the model from low energy obser- 
vations, let us discuss its embedding into the SU(2)l x SU(2)r x S77(4) c model [TO]. The 
version of the model relevant to our discussion is not the one in the original Pati-Salam 
paper but rather the one considered in Ref. [11] . In this model [11], symmetry breaking 
from SU{2)r x S77(4) c to U(l)y x SU(3) C is implemented by the Higgs fields belonging to 
the representation A R (1, 3, l"0) © A L (3, 1, fO) under the SU{2) L x SU(2) R x SU(4) e group. 
Decomposing this field under the standard model group SU(2)l x U(1)y x SU(3) c gives the 
various fields in the model: 

A R (1, 3, 10) = A u c u c(l, +|, 6*) © A uC(i c(l, +|, 6*) © A d c rfo (l, 6*) © A^l, +|, 3*) 

2 2 8 

©A d c i/C (l, --, 3*) © A u c e c(l, --, 3*) © A d c e c(l, --, 3*) 

©A^ c (l, 0, 1) © A e c„ c (l, -2, 1) © A eCe c(l, -4, 1) . (2) 
The Yukawa Lagrangian of this model is given by 

d = fi^fC-^T-K^f +(R^L)+ H^t^ + h.c. (3) 

I U\ U2 M3 ^\ 

where $ = . The / couplings generate Majorana neutrino masses, while 

the couplings denoted H a generate the Dirac masses for fermions. Comparing Eq. ([3]) with 
Eq. (0Q), we see that exactly the same interactions are present in both cases. The S field of 
Eq. ([T]) is the A u c u c whose vev breaks the gauge group of our model down to the SM gauge 
group [11]. We assume that the scale of this symmetry breaking is anywhere between 1 - 
100 TeV so that both the right-handed neutrinos as well as the gauge bosons belonging to 
SU(4) C /SU(3) C have masses around these values. 

We note that while breaking the gauge symmetry only by the A(l, 3, 10) makes the Wr 
mass scale, B — L breaking scale vrl, and the SU(4) C breaking scales all equal, and would 



also relate the Wr mass with the Z' mass. Some of the constraints we derive below require 
that Mw R 3> vbl- This can be achieved by including a (1,3, 1) Higgs field to break the 
symmetry, which will generate Wr mass, but not Z' mass and decouple Wr mass from vbl- 
Our results will be valid in the presence of such (1,3,1) Higgs fields, or in their absence. In 
the latter case, the common scale of B — L symmetry breaking will be required to be > 100 
TeV or so. 

The Higgs fields belonging to the A multiplet will have the following mass pattern: 
Arfc d c, A u c d c will have mass near a TeV, whereas A u c u c will have mass near 100 TeV. Such a 
mass pattern is consistent, since as noted above, we could have the SU (4) C /SU (3) c xU(1)b-l 
boson and the Wr mass different from the Z' mass and the right-handed neutrino mass scale. 

One important point to note is that due to SU (4) c gauge symmetry, all three couplings 
in Eq. (TjQ) become equal to each other i.e., hij = gij = ly = fij of Eq. (J3]). 

In general the neutrino mass in this model is given by a combination of type I and type 
II seesaw contributions: 



The coupling matrix / that appears in the neutrino mass formula above is related to the 
diquark couplings, which lead to FCNC effects, n — n oscillations as well as the baryon 
asymmetry. They are therefore very highly constrained. 

In what follows, we assume that M® aa,c = or very small by an appropriate choice of the 
Yukawa couplings of $ a ~ (2, 2, 1) and $ i5 ~ (2, 2, 15) fields (see Eq. ©). The details of this 
are not relevant to the main point of our paper. Since Mr, depends on the same parameters 
as the quark and charged lepton masses, it is useful to point out that setting Mb, = does 
not lead to any conflict with realistic fermion mass and mixing patterns. To see this note 
that the two bi-doublet fields $i and $15 are both complex scalars; therefore each field will 




(4) 
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have two independent Yukawa couplings with fermions: 

£ Yukawa = Y^l^r + Y^l^r + Y 15 1p L $ 15 l/j R + Yl 5 ^L$15^i? + h.C. (5) 

where ipL ~ (2, 1,4) and ipR ~ (1,2,4*) fermions, and ~ t 2 &*t 2 . From Eq. (jSJ), it follows 
that the Dirac mass matrices of the up quark, down quark, charged lepton, and the neutrino 
are all independent. 

Once we set M® nac = 0, we can directly link the neutrino mass matrix to the coupling 
matrix /. The advantage of this is that the requirement of adequate baryogenesis as well 
as consistency with FCNC and other constraints fix not only the neutrino mass matrix, 
but also the mass spectrum of the theory. The FCNC constraints come from the fact that 
A u c uC , A^c d c, A u c d c fields have masses in the multi-TeV range and can lead to sizable K° — 
K°, D° — D° and B® s — B® s mixings. They in turn severely constrain the pattern of the 
Yukawa couplings /y and thereby the neutrino mass matrix. 

3 Low energy constraints on the model 

In this section, we discuss the tree level flavor changing neutral current contributions to 
processes such as K — K, B^ 8 — B^ s ^ D — D mixings from the di-quark Higgs field exchanges. 
We have to make sure that they are not in conflict with observations. One cannot make the 
di-quark scalar masses very large to satisfy the FCNC constraints, since successful post- 
sphaleron baryogenesis requires the masses of at least two of these scalars to be not more 
than about a TeV. Similarly, the doubly charged scalar bosons from the same multiplet will 
contribute to rare processes such as fi — > 3e via tree level diagrams. Neutrino oscillation 
data, on the other hand, suggest specific form of the / matrix. We need to examine if 
these dual requirements can be simultaneously met. We have found that indeed this can be 
satisfied, but only for an inverted mass hierarchy spectrum for the neutrinos. 
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To discuss the constraints on the couplings fy and masses of A uV , A^c, A u c d c implied 
by these considerations, we first note that above the SU(4) C scale, all couplings to diquarks 
and dileptons are given by a single matrix /^-. The form of this matrix can be specified in any 
basis without loss of generality and we specify them in the basis in which the down quarks 
are mass eigenstates. In this basis, the fij couplings split up into the following depending 
on which quarks they couple to: f dd , fud and f uu , where f dd indicates the coupling to d c d c , 
etc. Assuming for simplicity that CP is not broken by the vacuum expectation values of the 
bidoublet fields (so that the left-handed and right-handed CKM matrices are equal to each 
other), we get 

fud = ^ CKM fdd 

fuu = UcKMfddUcKM 

fuu = U l f dd Uj = f ee (6) 

where Uqkm is the quark rotation matrix and U\ is the matrix that makes the charged 
leptons diagonal. Clearly, it is f vv which determines the neutrino mass matrix in the type 
II seesaw case. 

In this basis, first there are constraints from flavor changing processes such as K — K, 
B Sjd — B s d and D — D mixings. Below we list the constraints [13] and their implications for 
the parameters of the model: 

K°(ds) - K°(ds) : fddAlf ^ < 3.3 x 1(T 6 (7) 

B° s (sb) - B° s (sb) : fdd ' 22f ^ < 2.0 x 10" 4 (8) 
[m AdCdC /TeV] 2 

B° d (db) - B° d (db) : fddAlf % 3 " & 7.6 x 1(T 6 (9) 

[m AdCdC /TeV\ 2 

D°(uc) - DHuc) : ; fuuAlfu 2 ,2 L n £ 2 x 10" 6 (10) 



[m AuCuC /TeV] 
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In addition, lepton family number violating modes [H] such as [i — > 3e imply 

fee, 11 fee, 12 < q q -, n _ 5 ^\ 

[m A++ /TeV]2 ~ J d X 10 • (U) 

This can be satisfied by requiring the A ++ mass to be in the 100 TeV range for our choice 
of /i2,n as we see below. The constraints from the various r decay modes can then be easily 
satisfied for this limit on the A ++ mass and we do not give those constraints here. 

Another constraint on the parameters of the theory comes from the present limits on 
n — n oscillation period. r n _ fi > 10 8 sec. [T5l [T6] implies that the strength G n _ n of the 
AB = 2 transition is < 10~ 28 GeV~ 5 . In a generic model of this type, n — n oscillations arise 
from the tree diagram in Fig. [2] (see Sec. 5) and we find that 

„ ^1 (S) fi d ii/uu,ll A 2 (S) fdd,llfud,ll / ln -28 n \ T -h (m\ 

G n -f> — a + n a ■ — < 10 Gev . 12 

Mi Ml Ml Mi y ' 

We discuss this further in our model in Sec. 6. 

The nontrivial aspect of this model is that the same set of parameters responsible for 
baryogenesis is directly related to neutrino masses and mixings and must be such that they 
satisfy the strong FCNC constraints listed above. Note that we cannot suppress the FCNC 
effects by simply raising the masses of A^ crfc , A u c d c particles since in that case we cannot 
satisfy the desired constraints for adequate baryogenesis. 

4 Inverted neutrino mass hierarchy from the FCNC 
constraints 

In this section, we address the question of how we satisfy these constraints and yet obtain 
concordance with neutrino oscillation observations. It turns out that if we choose fdd matrix 
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as (in a basis where down quarks are mass eigenstates) 



/ 0.95 



dd 



\ 



1 

0.95 0.01 
V 1 0.01 -0.0627357 / 



(13) 



then for TeV scale A^c an d A uCrfc and 100 TeV mass for A u c u c, we can satisfy all the 
hadronic constraints. Such a choice will automatically satisfy K — K and Bd yS — Bd yS mixing 
constraints, owing to the zeros in the diagonal entries. In the leptonic sector, as already 
noted, the most stringent constraint comes from jj, — > 3e and it requires that A ++ mass also 
be of order 100 TeV or so. 

Then we use the following unitary transformation to rotate fdd to get the neutrino mass 
matrix: 



I cos 6 sin \ 
— sin cos 

V o 01/ 



(14) 



with = 0.23. This gives 

/ 0.421751 0.85125 0.975946 \ 

fu = UifddUf = 0.85125 -0.421751 -0.218241 (15) 

V 0.975946 -0.218241 -0.0627357 ) 

This matrix must be multiplied by the SU(2)l triplet vev vl (which is much smaller than 

the electroweak vev in Type II seesaw) to give the neutrino mass matrix. For vl = 0.03 eV, 

we get the three neutrino masses to be 

mi = 0.0410 eV, m 2 = -0.0417 eV, m 3 = -0.0012 eV 

which yields the correct solar-to-atmospheric neutrino mass square difference Am5; tm /Am;? olar : 
30 and the PMNS mixing angles 6 12 = 35.6°, 6 23 = 46° and 13 = 8°. This value of 13 is 
observable in the ongoing (Double CHOOZ) [17] and planned (Daya Bay) [18] experiments 
and should provide a test of the model. Secondly we predict that the neutrinoless double 
/?-decay experiments should observe a Majorana neutrino mass at the 10-20 meV level which 
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is perhaps within reach of the next round of neutrinoless double /3-decay experiments. In 
Fig. EH we present two scatter plots that display the preference of oscillation parameters in 
our model. We clearly see the lower bound on the #13 from them. 




0.015 0.02 0.025 0.03 



0.35 
0.34 
0.33 
0.32 
0.31 
0.3 
0.29 

0.015 0.02 0.025 0.03 0.015 0.02 0.025 0.03 

sin 2 8 !3 sin 2 e 13 sin 2 8 13 

Figure 1: We give the predictions for neutrino oscillation parameters for the allowed ranges 
of the di quark scalar couplings in our model. Note the lower limit on the #13 of about 0.1. 

We will see below that this form of the / matrix satisfies the baryon asymmetry con- 
straints as well as the n — n constraints. 

5 Origin of matter 

Before proceeding to the discussion of how baryon asymmetry arises in this model, let us first 
sketch the cosmological sequence of events starting at the SU (4) c scale that leads up to this. 
For temperatures above the 577(4) c scale of about 100 TeV, there is no B — L violation. The 
sphalerons are active and therefore erase any pre-existing B + L asymmetry in the universe. 
So if there was a primordial GUT scale generated baryon asymmetry that conserved B — L 
(like that in most SU(5) and some SO( 10) models), it will be erased by sphalerons. Any 
baryon asymmetry residing in B — L violating interactions will however survive. 

Below the S77(4) c scale, B — L violating interactions arise e.g. SS — > e + e", and will 
be in equilibrium together with the AB = 2 interactions. So together they will erase any 
pre-existing baryon or lepton asymmetry. Thus in models of this kind, baryon asymmetry 
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of the universe must be generated fresh below the sphaleron decoupling temperature. 

In order to sketch how fresh baryon asymmetry arises in our model, we assume the 
following mass hierarchy between the S r field and the A^c, A u c u c, A u c d c fields: 

m t < M s (~ 500 GeV) < M AdCdC ~ M A „ ede (~ 1 TeV) < M AuCuC (~ 100 TeV), 

where m t is the top quark mass. 

Between 1 < T < 100 TeV, the AB = 2 interaction rates go like 

T(AB = 2) ~ J^T (16) 

and are therefore in equilibrium if some of the /^'s are above 0.3 as in our case. 

Below T ~ 1 TeV, the AB = 2 processes such as the decay S r — > 6q c +6q c , (q c , q c )+ S r — > 
5(g c , q c ) occur at a rate given by 

Y(AB = 2) ~ 100/ ^ 12 (17) 
1 ; (2tt) 9 (6M) 12 1 ; 

where M ~ TeV, the average mass of the A^^, A u c d c particles which are still in equilibrium. 
The A u c u c is about 100 TeV and hence its contribution to these processes is more suppressed 
compared to that of A d c d c, A u c d c. This decay then goes out of equilibrium somewhat below 
the TeV temperature range. One impact of this is that these interactions being in equilibrium 
above T ~ TeV erase any pre-existing baryon asymmetry as discussed above. 

By the time the universe cools to temperature near or slightly below Mg, its decay 
channels can start if the rates are faster compared to the Hubble expansion rate. Let us 
therefore estimate the various decay rates: 

There are four decay modes which are competitive with each other: (i) S r — > 6g c ; (ii) 
S r — > Zf c f c ; (iii) S r — > ZZ and (iv) S r — > re[j] We discuss them below. 



1 The S r — » W + W is suppressed by Wl — Wr mixing parameter which can be adjusted to be small. 
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d c (ti c ) 



Figure 2: Tree level diagrams contributing to S r decays into 6 anti-quarks. There are other 
diagrams where S r decays into 6 quarks, obtained from the above by reversing the arrows of 
the quark fields. 

(i) S r — > 6q c : The diagram for this is given in Fig. [2J Since Ms 3> m t , in its decay all 
modes will participate. Including all the modes, we find the decay rate to be: 

T(s 6<f) ~ 36 Wm^l (1 8) 

(b r -»• bq ) _ (2?r)9 (6A ^ )12 (18) 

where we have chosen Ai = A2 = A ~ 0.1. Taking as an example a typical set of parameters 
Ma — 2Ms ~ 1 TeV and taking the parameters for the / matrix elements from Eq. ( |T3l) . we 
get T(S r -> Qq c ) ~ 7.5 x 10' 17 GeV. 

(ii) S r — > Z + f c f c : This arises from the SS coupling to Z'Z' with one of the Z's mixing 

with Z' (Fig. 3) and the virtual Z' decaying to f c f c . This occurs only for T < v w k- This is 

because for T > v w k, Z — Z' mixing disappears. Below the electroweak symmetry breaking 

temperature, this mixing denoted below by gzz> becomes effective and is given by 

g 2 cos 2 9 w (M z \ 2 
Qzz> = , TT~ v B l (19) 



Vcos2^ \M Z 

which leads to the new S r decay mode (Fig. 3) given above. This decay rate is given by 

7.0 x 10- 2 GeV 2 



T(S r -> Zffc 



M S M% 



M A M; - Ml (6M S 4 - 19M 2 M| + 28M. 



, M s + \/M 2 — M 2 

-3M|(M 2 + 4Mf) log 



(20) 



13 



Figure 3: Feynman diagram for S r decay to Zf c f c vis Z — Z' mixing 

For our choice of parameters and Mz> ~ 100 TeV, we find that T(S r — > Zf c f c ) ~ 6.6 x 10~ 18 
GeV and is therefore slower than the 6q° decay rate. 

(iii) S r — > ZZ : This decay mode arises from Z — Z' mixing with the decay width given 

by 

T(S ZZ) - 9hMl (l AM2 A 
where the S r ZZ vertex is given by 

9zz = ^g 2 cos 2 9 w v B L {jf^j ■ ( 22 ) 

For M Z ' ~ 100 TeV and M s = 500 GeV, we get F(S r -> ZZ) ~ 2.2 x 1(T 19 GeV which is 
much smaller than the Qq c decay rate. 

In Fig. HI these decay rates are plotted against the mass of the scalar field for various 
values of vbl- Note that Mz> is related to Vbl as follows: 

M 2 ; _ 2 9 2 V 2 BL COS 2 W 

z ' cos 26 w 



1 - + Z 



Ml 



Ml 



(21) 



(iv) S r — > r + e: This decay mode arises from the Feynman diagram in Fig. and its 
rate can be estimated to be: 

1 ^ r - r + ej _ 127r(167r2)232M ^ • l^J 
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10 " ^ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — " 

200 400 600 800 1000 

M s (GeV) 

Figure 4: The S r — > Zf c f c (thin solid lines) and S r — ► ZZ (thin dashed lines) decay rates 
for various values of vbl (in TeV). The thick solid lines correspond to the S r — > 6q c decay 
rates for two typical values of r = M& uCdC dCdC /M s . We see that for vbl > 40 TeV, the six 
quark decay mode dominates for a large range of Ms- 




Figure 5: S r — > er decay 



This width is estimated to be T(S r —>■ r + e) ~ 9 x 10 20 GeV. Therefore this is also much 
smaller than the decay rate to six quark modes. 

At the time the universe has a temperature of ~ Ms or slightly below so that it is out of 
equilibrium from the cosmic soup, the Hubble expansion rate is ~ ^JglM'gjMpi ~ 2.5 x 10~ 12 
GeV implying that all the above decay modes are out of equilibrium. Since the decay 
rate remains constant below T ~ M$, but the expansion rate of the universe is slowing 
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down as it expands, there will come a time (or temperature Td) when the dominant decay 
T(S r — > 6q) ~ H{Td). At that point the S r particle will start decaying and produce the 
baryon asymmetry as in Ref. [8]. 

At this temperature which is far below the masses of the A uC(JC , A^c particles, the decay 
processes A q c q c — > q c q c being very fast have depleted all the diquark Higgses and have left 
only the S r particles to survive along with the usual standard model particles. The primary 
decay modes of S r are S r — ► u c d c d c u c d c d c and S r — > u c d c d c u c d c d c as already noted (Fig. EJ). 
Other decay modes are negligible as discussed. 

We have to make sure that the decay of S r starts below the sphaleron decoupling tem- 
perature and above the QCD phase transition temperature. To check if this indeed happens 
in our model, let us calculate the T^: 

rp „ [ 36A 2 (Tr[/t / ])3 M p 1 Ml 3 l 1/2 
d [ (2vr) 9 1.66^ 1/2 (6M A ) 12 . 

= 6.1 GeW^ (24) 

For M s ~ 500 GeV and M A ~ M AuCdC ~ M AdCdC ~ 1 TeV, we get T d ~ 2 GeV which is 
comfortably above the QCD phase transition temperature. 

It is worth emphasizing that if we increased the sextet scalar masses arbitrarily to satisfy 
the FCNC constraints, this will lower the Td to undesirable values below the QCD temper- 
ature. One may think that we could simultaneously increase the value of Ms but as we will 
see below, the magnitude of the baryon asymmetry goes inversely like the square of Ms and 
increasing it above 500-600 GeV will suppress the baryon asymmetry to a level below the 
observations. 

The calculation of the baryon asymmetry is same as in Ref. [8] and we do not repeat 
it here except to give a brief summary for the particular / texture in our model. First 
we note that since only /31 and /21 are the dominant contributions, there is a flavor factor 
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Figure 6: One loop vertex correction diagram for the B- violating decay S r — > 6q c . There 
are also wave function corrections involving the exchange of W ± gauge bosons, which are 
somewhat smaller. 

of 64 = (Tr [/'/] ) 3 i n the absolute decay width. Next we calculate the baryon asymmetry 
from the vertex correction via the W boson exchange (Fig. 6), which dominates the baryon 
asymmetry. To do this we note that the dominant contribution comes from making the / 
matrix complex (e.g. ^33 with a maximal complex phase and all other parameters real). This 
of course does not affect the neutrino fit discussed earlier. In terms of the flavor combinations 
that give the dominant contributions, they come from products like /fi/33 of which there are 
six combinations. This gives 



We have also assumed that Ms 3> m t . Note that if we increased Ms above 500 GeV or so, 
the generated baryon asymmetry will fall short of the observed values. 

This gives for the baryon asymmetry at T = T^. ee ~ (2 — 3) x 10~ 8 . To compare it with 
the observed t]b, we divide this by g*(200 MeV)/gf*(l eV) ~ 62.75/5.5 = 11.4 and apply 
an additional dilution factor of 0.25 (see discussion below) which gives us the desired value. 
Note that the observed value of the asymmetry is rj^ MB ~ 6 x 10 -10 [19J. 

Since the SV-particle decays far below its mass to generate the baryon asymmetry, we 



.vertex 
B 



Br 



a 2 6 Im [f^mtVtbmbf^mtVtbmb} 
4 (Tr[/t / ])3 M ^M| 



(25) 



17 



have to take into account the effect of its decay, which as we explain below amounts to 
a dilution of the original baryon asymmetry calculated. In order to estimate the dilution 
factor, we note that SV-decay will release all the energy in its mass to lighter relativistic 
particles which will thermalize with the rest of the cosmic fluid and in the process raise its 
temperature which will increase the entropy and hence dilute the net baryon asymmetry. 
Suppose the decay temperature is T d . Energy conservation then gives 



PS + Prc\\ Td — Prel\ T> , 

h* T *M s + —g m T* =-. 
vr 2 d 6 3(T d 30' 



or, ~^TlM s + -gj* = -g*T* . (26) 



Solving this one finds for the dilution factor d that 

d = H ^ °- 32 ^ T > (27 ) 

T| 0.1271^ + 0.32^' 1 ; 

For M s ~ 500 GeV and T d ~ 1 GeV, d ~ 0.25. 

There is another factor coming from the fact that the baryon asymmetry generated is at 
T d ~ 1 GeV where g* = 62.75 whereas the measured value is at T rec where g* = 5.5. All 
these dilution factors have been taken into account in our estimate of final baryon asymmetry, 
which is in agreement with observations. 

6 Prediction of observable neutron-antineutron oscil- 
lation time 

In this section, we discuss the prediction of the model for neutron-anti-neutron oscillation. 
In order to estimate the n — n oscillation, let us first recall that the only contribution to 
this process comes from the right-handed sector since the vev of the A u c u c is in the 100 TeV 
range and that of its left-handed counterpart is in the eV range. 
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There are two types of contributions to n — n oscillation from the right-handed sector: 
the one involving two d c d c type and one u c u c type bosons of the right-handed sector and 
another which involves two u c d c and one d c d c type A boson. Since the diagonal 11 and 22 
entries of fdd are close to zero, the first contribution is actually much smaller than the second 
one. The second type generates an effective operator of the form u c d c b c u c d c b c . To get n — n 
oscillation, we will have to change the two b c quarks to two d c quarks by second order weak 
interactions (see Fig. [7]). 





" L 



Figure 7: Loop diagram for n — n oscillation 



From Fig. [71 we see that the six quark operator for n — n oscillation has 7^'s in the 
Lorentz structure i.e. (uR^^di ) 2 d c d c . This is different form for the AB = 2 operator whose 
matrix element has not been evaluated in the literature. For the strength of this operator 
we estimate 



Gn- 



fud,nfud,i3fdd,i3^VBL g\ 2 d m 2 b m* 

-lUg( 

itt c dC — t± d c d c v-~- / ■■ b w m w 



(28) 



This gives G n _ n ~ 10~ 30 GeV -5 . Taking the hadronic "dressing" of quark to hadrons to be a 
factor of 10~ 4 , we estimate n — n transition time of 10 9-10 sec. given the uncertainties in the 
parameters. The present lower limit on this transition time is 10 8 sec. from the Grenoble 
experiment [15J as well as from nuclear decay experiments [IB] . Our predicted value is 
accessible to current experiments under discussion at DUSEL as well as other facilities [20J . 
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7 Other implications and tests of the model 



(i) As noted in Sec. 2, a crucial prediction of our quark-lepton unified model of post-sphaleron 
baryogenesis is that neutrinos must be Majorana fermions and exhibit an inverted mass 
hierarchy form with large value for # 13 . This should be testable in long base line experiments 
as well as the ongoing and planned reactor experiments searching for #13 and neutrinoless 
double beta decay searches [21]. It is perhaps worth noting that there is indication of a 
non-zero #13 from already existing neutrino oscillation data [22J. 

(ii) Our theory is also testable in collider experiments such as the LHC since we have colored 
diquark scalar fields with masses in the TeV range. It is clear from the form of the f ud matrix 
that in a pp collision, the valence quarks in the two protons could produce the A u c d c field 
which could then decay to t+ jets. This could either be an s-channel single production [23] 
or Drell-Yan pair production [24] . The s-channel process will have a resonant enhancement 
which can give a signal above the standard model background. The Drell-Yan pair production 
could give signals of type bbl ± l ± jj+ missing Et- Unlike the s-channel process, the Drell- 
Yan pair production has the advantage of not being dependent on the specific flavor texture 
of the two quark couplings / and is promising for color sextet masses upto a TeV [24J. 
It would therefore be important to search for these particles at LHC. Their discovery will 
signal a completely different direction for unification beyond the standard model than the 
conventional SUSY-GUT theories. 

(iii) In our model, since there is a mass hierarchy between the A uC(JC , and A u c u c masses 
i.e. M& uCdC , M^ dCdC <C Ma u c„ c , a one loop level box graph induced by the trilinear coupling 
\vBLA. u cd,cA u c d c&. d c d c will induce a quartic coupling (A* C(JC A uC d c ) 2 coupling with a strength 
\ e ff ~ — ^ a/^ bl J ■ This can lead to color breaking unless A e // < 1 — 2. This can be 
satisfied by lowering the vbl scale to about 50 TeV with A ~ 0.1. In order to reconcile this 
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lower value with constraints from \i — > 3e, we can introduce a multiplet of type (1, 3, 1) with 
a vev in the 100 TeV range which gives mass to the Wr and the A ++ . This vev decouples 
the B — L breaking scale vbl from the masses of the Wr and A ++ fields. In this case, one 
can keep the B — L breaking scale near 50 TeV while keeping the Wr and A ++ mass around 
100 TeV as required by the fi — > 3e and S r — > r + e constraints. 

(iv) Finally, note that a priori in the model there could be a coupling of type A^ATr[$^$], 
which will induce a S r decay to two SM Higgs fields. We assume that this parameter is very 
small . This assumption could be justified in supersymmetric extensions of the model where 
such terms are forbidden by holomorphy of the superpotential. 

8 Conclusion 

In summary, we have pointed out that the post-sphaleron baryogenesis mechanism proposed 
in Ref. 0,19] can be naturally embedded into a 100 TeV scale quark-lepton unified SU{2)i x 
SU(2)r x SU(4) C model. If we further assume that the neutrino masses in this model arise via 
the type II seesaw mechanism, then the couplings responsible for baryogenesis and neutrino 
masses get intimately linked to one another. In this case adequate baryogenesis predicts that 
neutrino mass ordering must be inverted with large #13, a prediction that can be tested in 
ongoing neutrino experiments. 
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04ER41306, and that of RNM and BD is supported by the National Science Foundation 
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